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^ ' The perturbative QCD approach is based on fcy factorization, including the Sudakov 

pH ' form factors so that to avoid the endpoint singularity. In this approach, we calculate the 

charmless B decays like B — * tttt decays etc. to produce the right number of branching 
ratios and also CP asymmetry parameters. For final states with at least one charmed 
meson, like B — > Diz decays, our results also agree with the experiments. 
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The hadronic B decays are important for the CKM matrix elements measurements 

and CP violation detection. Understanding nonleptonic B meson decays is crucial 

for testing the standard model, and also for uncovering the signal of new physics. 

, The simplest case is two-body hadronic B meson decays, for which Bauer, Stech and 

?H ' Wirbel (BSW) proposed the naive factorization assumption (FA) in their pioneering 

rrt 1 rn |Ai<>| 

workW. Considerable progress, including generalized FA - and QCD-improved FA 

(QCDF) SI, has been made since this proposal. On the other hand, technique to 

analyze hard exclusive hadronic scattering was developed by Brodsky and Lepage ^ 

based on collinear factorization theorem in perturbative QCD (PQCD). A modified 

framework based on kx factori zation theorem was then given in 1^13 and extended 

to exclusive B meson decays fa EMM The infrared finiteness and gauge invariance 

of kr factorization theorem was shown explicitly in El Using this so-call ed PQCD 

11 911 4.I 

approach, we have investigated dynamics of nonleptonic B meson decays 1 1 1 I 
Although the predictions of branching ratios agree well with experiments in most 
cases, there are still some theoretical points unclear in FA and QCDF. First, it relies 
strongly on the form factors, which cannot be calculated by FA itself. Secondly, the 
generalized FA shows that the non-factorizable contributions are imp ortant in a 
group of channels, which can not be done reliably in FA and QCDF 12121 The reason 
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of this large non-factorizable contribution needs more theoretical studies. Thirdly, 
the strong phase, which is important for the CP violation prediction, is quite sensi- 
tive to various approaches. The mechanism of this strong phase is quite different for 
different method, and give quite different results. The recent experimental results 
can make a test for the validity of these approaches. 

2. Formalism of PQCD Approach 

In this section, we will introduce the idea of PQCD approach. The three scale PQCD 
factorization theorem has been developed for non-leptonic heavy meson decays 
based on the formalism by Brodsky and Lepage^, and Botts and Sterman^. In the 
non-leptonic two body B decays, the B meson is heavy, sitting at rest. It decays 
into two light mesons with large momenta. Therefore the light mesons are moving 
very fast in the rest frame of B meson. In this case, the short distance hard dynamic 
dominates the decay amplitude. The reasons can be ordered as: first, because there 
are not many resonance near the energy region of B mass, so it is reasonable to 
assume that final state interaction is not important in two-body B decays. Second, 
With the final light mesons moving very fast, there must be a hard gluon to kick 
the light spectator quark (with small momentum) in the B meson to form a fast 
moving light meson. So the dominant diagram in this theoretical picture is that one 
hard gluon from the spectator quark connecting with the other quarks in the four 
quark operator of the weak interaction. Unlike the usual factorization approach, the 
hard part of the PQCD approach consists of six quarks rather than four. We thus 
call it six-quark operators or six-quark effective theory. There are also soft (soft 
and collinear) gluon exchanges between quarks. Summing over those leading soft 
contributions gives a Sudakov form factor, which suppresses the soft contribution 
to be dominant. Therefore, it makes the PQCD reliable in calculating the non- 
leptonic decays. With the Sudakov resummation, we can include the leading double 
logarithms for all loop diagrams, in association with the soft contribution. 

There are three different scales in the B meson non-leptonic decay. The QCD 
corrections to t he f our quark operators are usually summed by the renormalization 
group equation ^^-1 This has already been done to the leading logarithm and next- 
to-leading order for years. Since the 6 quark decay scale mj, is much smaller than 
the electroweak scale mw, the QCD corrections are non- negligible. The third scale 
1/6 involved in the B meson exclusive decays is usually called the factorization 
scale, with b the conjugate variable of parton transverse momenta. The dynam- 
ics below 1/6 scale is regarded as being completely non-perturbative, and can be 
parameterized into meson wave functions. The meson wave functions are not cal- 
culable in PQCD. But they are universal, channel independent. We can determine 
them from experiments, and it is constrained by QCD sum rules and Lattice QCD 
calculations. Above the scale 1/6, the physics is channel dependent. We can use 
perturbation theory to calculate channel by channel. 

Besides the hard gluon exchange with the spectator quark, the soft gluon ex- 



I 



I 



1, 2008 9:51 WSPC/INSTRUCTION FILE lucd-guangzhou 
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changes between quark lines give out the double logarithms In 2 (Pb) from the overlap 
of collinear and soft divergence, P being the dominant light-cone component of a 
meson momentum. The resummation of these double logarithms leads to a Sudakov 
form factor exp[— s(P, b)], which suppresses the long distance contributions in the 
large b region, and vanishes as b > 1/Aqcd- So this term includes the double in- 
frared logarithms. The corresponding Sudakov factor can be derived in PQCD as 
a function of the transverse separation b and of the momentum fraction x carried 
by the spectator quark 121 . The Sudakov factor suppresses the large b region, where 
the quark and antiquark are separated by a large transverse distance and the color 
shielding is not so effective. It also suppresses the x ~ 1 region, where a quark car- 
ries all of the meson momentum, and inten ds to emit real gluons in hard scattering. 
The Sudakov factors from kx resummation for the B and Z?W mesons are only 
associated with the light spectator quarks, since the double logarithms ar ise from 
the overlap of the soft and mass (collinear) divergences. It is shown in refP^ that 
e~ s falls off quickly in the large b, or long-distance, region, giving so-called Sudakov 
suppression. This Sudakov factor practically makes PQCD appro ach applicable. For 
the detailed derivation of the Sudakov form factors, see refl ^^ l 

With all the large logarithms resummed, the remaining finite contributions are 
absorbed into a perturbative b quark decay subamplitude H(t). Therefore the three 
scale factorization formula is given by the typical expression, 



C(t) x H(t) x ®(x) x exp 



-s(P,b)-2 f % q (a s {fi)) 
Ji/b A 1 



(1) 



where C(t) are the corresponding Wilson coefficients, <j>(x) are the meson wave 
functions and the variable t denotes the largest mass scale of hard process H, that 
is, six-quark effective theory. The quark anomalous dimension ^ q — —a s /Tt describes 
the evolution from scale t to 1/6. Since logarithm corrections have been summed by 
renormalization group equations, the above factorization formula does not depend 
on the renormalization scale /i explicitly. 

As shown above, in the PQCD approach, we keep the kx dependence of the 
wave function. In fact, the approximation of neglecting the transverse momentum 
can only be done at the non-endpoint region, since kx <C k + is qualified at that 
region. At the endpoint, k + — > 0, kx is not small any longer, neglecting kx is a very 
bad approximation. By, keeping the kx dependence, there is no endpoint divergence 
as occurred in the QCD factorization approach, while the numerical result does 
not change at other region. Furthermore, the Sudakov form factors suppress the 
endpoint region of the wave functions. Recently another type of resummation has 
been observed. The loo p correction to the weak decay vertex produces the double 
logarithms a s In 2 xn ^1 Using the wave functions from light-cone sum rules, at the 
endpoint region, these large logarithms are important, they must be resummed. 
The threshold resummation for the jet function results in Sudakov suppression, 
which decreases the contribution of endpoint region of wave functions. Therefore, 
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Fig. 1. Diagrams for B — > P1P2 decay in perturbative QCD approach. The factorizable diagrams 
(a),(b), non-factorizable (c), (d), factorizable annihilation diagrams (e),(f) and non- factorizable 
annihilation diagrams (g),(h). 



the main contributions to the decay amplitude in PQCD approach comes not from 
the endpoint region. The perturbative QCD is applied safely. 

The main input parameters in PQCD are the meson wave functions. It is not a 
surprise that the final results are sensitive to the meson wave functions. Fortunately, 
there are many channels involve the same meson, and the meson wave functions 
should be process independent. In all the calculations of PQCD approach, we follow 
the rule, and we find that they can explain most of the measured branching ratios 
of B decays. For example: B — > irir decays, B — > irp, B — > iru> decays B — > Kir 
decays ^21 b — > KK decays the form factor calculations of B — > n, B — > 
B -> Krft) decays^, B — > Kcj) decays^ etc. 

We emphasize that nonfactorizable and annihilation diagrams are indeed sub- 
leading in the PQCD formalism as Mb — > 00. This can be easily observed from the 
hard functions in appendices of ref.™. When Mb increases, the B meson wave 
function enhances contributions to factorizable diagrams. However, annihilation am- 
plitudes, being independent of B meson wave function, are relatively insensitive to 
the variation of Mg. Hence, factorizable contributions become dominant and an- 
nihilation contributions are subleading in the Mb — > 00 limit ^2. Although the 
non-factorizable and annihilation diagrams are subleading for the branching ratio 
in color enhanced decays, they provide the main source of strong phase, by inner 
quark or gluon on mass shell. This mechanism of strong phase is negligible in the 
PQCD approach, since it is at next-to-leading order 0(a s ) corrections. In fact, the 
factorizable annihil ation d iagrams are Chirally enhanced. They are not negligible 
in PQCD approach 

In the PQCD framework based on fcy factorization theorem, an amplitude is 
expressed as the convolution of hard b quark decay kernels with meson wave func- 
tions in both the longitudinal momentum fractions and the transverse momenta of 
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QCD in hadronic B decays 5 



Table 1. Predicted B — > Dir decay amplitudes in units of 10 2 GeV, and branching 
ratios in units of 10 — 3 . 



Quantities 


C D = 0.6 


C D = 0.8 


C D = 1.0 


Data 


Mi 


6.39 - 0.35i 


6.88 - 0.38i 


7.35 - 0.40? 




M 2 


-1.53 + 1.48* 


-1.49 + 1.48* 


-1.45 + 1.45? 




M 3 


8.56 - 2.45i 


8.99 - 2.47i 


9.40 - 2.46i 




B(B U -> D+7T") 


2.37 


2.74 


3.13 


3.0 + 0.4 


B(B° -> D°tt°) 


0.26 


0.25 


0.24 


0.29 + 0.05 


B(B~ -» D°tt-) 


4.96 


5.43 


5.91 


5.3 + 0.5 



partons. In the B — > Dtt like decays with at least one heavy meson in the final 
states, our PQCD formulas are derived up to leading-order in a s , to leading power 
in rri]j/mB and in A/ron, and to leading double- logarithm resummations. 

3. Numerical Results and Discussion 

The PQCD predictions for each term of the B — > Dtt decay amplitudes are 
exhibited in Table [T] The theoretical uncertainty comes only from the variation of 
the shape parameter for the D meson distribution amplitude, 0.6 < Cd < 1.0. It is 
expected that the color-allowed factorizable amplitude / w £ext dominates, and that 
the color-suppressed factorizable contribution fo^iat is smaller due to the Wilson 
coefficient C% + CijN c ~ 0. The color-allowed non-factorizable amplitude .M e xt is 
negligible: since the pion distribution amplitude is symmetric under the exchange 
of X3 and I — X3, the contributions from the two diagrams Figs.QJc) andfljd) cancel 
each other in the dominant region with small x-i . It is also down by the small Wilson 
coefficient C%/N c . For the color-suppressed non-factorizable contribution .Mint, the 
above cancellation does not exist in the dominant region with small X3, because 
the D meson distribution amplitude cpuix^) is not symmetric. Furthermore, .Mint, 
proportional to C^/iVc ~ 0.3, is not down by the Wilson coefficient. It is indeed 
comparable to the color-allowed factorizable amplitude /^cxt, and produces a large 
strong phase. Both the factorizable and non-factorizable annihilation contributions 
are vanishingly small. 

The predicted branching ratios in Table Q] are in agreement with the averaged 
experimental data^l We extract the parameters a± and a-i from our calculations of 
amplitudes. That is, our a± and ai do not only contain the non-factorizable ampli- 
tudes as in gene ralized FA, but the small annihilation amplitudes, which was first 
discussed in^. We obtain the ratio |a 2 /ai| ~ 0.43 with 10% uncertainty and the 
phase of relative to a\ about Argia^jai) ~ —42°. Note that the experimental 
data do not fix the sign of the relative phases. The PQCD calculation indicates 
that Arg(aij a\) should be located in the fourth quadrant. It is evident that the 
short-distance strong phase from the color-suppressed nonfactorizable amplitude is 
already sufficient to account for the isospin triangle formed by the B — > Dir modes. 
Hence, it is more reasonable to claim that the data just imp ly a large strong phase, 
but do not tell what mechanism generates this phase ^4 From the viewpoint of 
PQCD, this strong phase is of short distance, and produced from the non-pinched 
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Table 2. Direct CP asymmetries calculated in FA, QCDF and PQCD 
for B — > ■kit and B — > Ktt decays together with the experimental 
results at percentage. 



Quantities 


FA 


QCDF 


PQCD 


Data 


B" -> tt+tt- 


-5 ±3 


-6± 12 


+30 ± 20 


+37 ± 11 


B° -> w+K- 


+ 10 ±3 


+5 ±9 


-17 + 5 


-10.9 + 1.9 


B + -» K°w+ 


+1.7±0.1 


+1±1 


-1.0 + 0.5 


-2.0 + 3.4 




+8 ±2 


+7 ±9 


-13 + 4 


+4 + 4 



singularity of the hard kernel. Certainly, under the current experimental and the- 
oretical uncertainties, there is still room for long-distance phases from final-state 
interaction. Other decays with one D meson in final states, like B — > D^K W 123 
B — > D^rf-'} and B — > D s (f>&^ etc. also agree with experiments. 

As discussed in the previous section, the strong phase generated from PQCD 
approach is quite different from the FA and QCDF approaches. The direct CP asym- 
metry is proportional to the sine of the strong phase difference of two amplitudes. 
Therefore the direct CP asymmetry will be different if strong phase is different. The 
predicted CP asymmetry by the three methods are shown in table[2] It is easy to see 
that the FA^and QCDFHTresults are quite close to each other, since the mechanism 
of strong phase is the same for them. Recently the two B factories measure at least 
two channels with non-zero direct CP asymmetry: B° — > ir + ir~ and B° — > ir + K~ 

which are shown in table [H It is easy to see that our PQCD results of direct 
CP asymmetry agree with the experiments. Although FA and QCDF are not 
yet ruled out by experiments, but the experiments at least tell us that the domi- 
nant strong phase should come from the mechanism of PQCD not QCDF. Charm 
quark loop mechanism, which gives the central value of strong phase in QCDF, is 
argued next-to leading order in PQCD. This argument is now proved by B factories 
experiments. 

4. Summary 

In the PQCD approach, the form factors are calculable, which are dominant by 
short distance contribution. By including the kx dependence and Sudakov sup- 
pression, there is no endpoint divergence. In the PQCD formalism non-factorizable 
amplitudes are of the same order as factorizable ones in powers of I/Mb, which 
are both 0(1/(MbAqcd))- The smaller magnitude of nonfactorizable amplitudes 
in color enhanced decays are due to the cancellation of the two non-factorizable 
diagrams. From the viewpoint of power counting, they are of the same order. In 
case of B — > Dtt decays, the cancellation is absent. The power counting changes, so 
that we can also calculate the non-factorizable contribution dominant processes. In 
this case, we give the right branching ratios for B — » D°ir° decay. 

The strong phase comes mainly from the annihilation and non-factorizable di- 
agrams in PQCD approach, which is quite different from the FA and QCDF ap- 
proaches. The experimentally measured direct CP asymmetry implies that PQCD 



February 1, 2008 9:51 WSPC/INSTRUCTION FILE lucd-guangzhou 



QCD in hadronic B decays 7 

gives at least the dominant strong phase than other approaches. This will be further 
tested by experiments. 
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